Magnetically induced chessboard pattern in the conductance of a Kondo quantum dot 
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We quantitatively describe the main features of the magnetically induced conductance modula- 
tion of a Kondo quantum dot - or chessboard pattern - in terms of a constant-interaction double 
quantum dot model. We show that the analogy with a double dot holds down to remarkably low 
magnetic fields. The analysis is extended by full 3D spin density functional calculations. Introduc- 
ing an effective Kondo coupling parameter, the chessboard pattern is self-consistently computed as 
a function of magnetic field and electron number, which enables us to quantitatively explain our 
experimental data. 



A quantum dot [1] with a finite net electron spin 
strongly coupled to its leads, enabling higher-order co- 
tunneling processes, is able to exhibit the Kondo effect 
[2,3]. The Kondo effect in quantum dots manifests it- 
self as an enhanced conductance in the Coulomb block- 
ade regime and occurs for temperatures and source-drain 
voltages below an energy scale set by the Kondo temper- 
ature. The first experimental results on the Kondo effect 
in quantum dots [4-6] were described with help of the 
spin- 1/2 Anderson impurity model. Assuming continu- 
ous filling of spin-degenerate single-particle levels, the dot 
is expected to either have total electron spin S = (for 
electron number N is even) or S = 1/2 (for TV is odd). 
The Kondo effect is therefore only expected for odd TV, 
hence giving rise to an "even-odd-effect" in the Coulomb 
valley conductance. A wide range of experiments, how- 
ever, has shown a clear deviation from this picture [7-16]. 
Particularly striking is the observation of a "chessboard 
pattern" in the dot conductance as a function of magnetic 
field, and gate voltage, V g [11,12,14-16]. Character- 
istic for this pattern is the alternation of high and low 
valley conductance regions as a function of B within the 
same Coulomb valley, i.e. for constant TV. In addition, 
the conductance also alternates when TV is changed by 
sweeping V g at constant B. The distinct regions in the 
V gi B plane of either high or low conductance are asso- 
ciated with the fields of a chessboard due to the similar 
appearance when the conductance is plotted in color scale 
(see Fig. lb). It has been experimentally shown that the 
enhanced conductance in certain Coulomb blockade re- 
gions can be ascribed to the Kondo effect, for both TV 
odd and even [11,12,14-16]. 

In this paper, we start by presenting our experimental 
data on a single lateral quantum dot clearly exhibiting 
the chessboard pattern. Next, we calculate the skele- 
ton of the chessboard pattern with a constant inter- 



action (CI) model of two capacitively interacting dots, 
formed by the two lowest Landau levels (LLs) - an anal- 
ogy (partly qualitatively) applied before at high magnetic 
fields [12,18-20] - and obtain the characteristic hexago- 
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FIG. 1. (a) Scanning electron micrograph of the device. 
The metal gates are yellow and the dot is indicated by a 
red circle. The ungated 2DEG has a mobility of 2.3 x 10 6 
cm 2 /(Vs) and an electron density of 1.9 x 10 15 m -2 at 4.2 
K. The nominal dot size is 320 x 320 nm 2 . (b) Color scale 
plot of the experimental linear conductance G through the dot 
as a function of B and V g i at 10 mK. The dotted hexagons 
highlight the shape of a few chessboard fields. The ellipses in- 
dicate some regions where suppression of the Coulomb peak 
occurs, (c) Calculated self-consistent potential landscape of 
the device. The white lines denote the contours of the metal 
gates. 
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nal, double dot (DD) stability diagram [17]. Interest- 
ingly, in this work we find that, due to small N and the 
shallow potential, higher LLs only start to fill below a few 
tenths of a tesla. Thus the DD analogy applies down to 
remarkably low magnetic fields. We finally perform full, 
3D spin density functional (SDF) calculations for our de- 
vice and introduce an effective Kondo coupling, derived 
entirely from the self-consistent results. Hence we can 
simultaneously calculate the electronic states of the dot 
and an estimate of the Kondo coupling, exhibiting the 
chessboard structure, and quantitatively explain some of 
its subtler features. 

Our quantum dot is shown in Fig. la. Metal gates 
are deposited on top of a GaAs/AlGaAs heterostructure 
with a two-dimensional electron gas (2DEG) 100 nm be- 
low the surface [17]. By depleting the 2DEG below the 
gates, the quantum dot is defined. Current can flow from 
the source (S) to the drain (D) contact. The electron 
number is varied by sweeping the left gate voltage, V g \. 
The SDF calculations show that our dot typically con- 
tains ~20 — 40 electrons. 

Figure lb shows a color scale plot of the linear con- 
ductance G through the dot as a function of B and V g \. 
Red (blue) corresponds to large (small) G (see scale in 
Fig. lb). For the most negative values of V g i the cou- 
pling of the dot to the leads is weak. This results in rel- 
atively sharp Coulomb peaks (red lines) and low valley 
conductance (dark blue regions). However, if V g i is in- 
creased, the valley conductance reaches considerable val- 
ues (~e 2 /h) in certain regions of the (B, V g i) plane. Most 
strikingly, the regions of low and high valley conductance 
alternate both along the V g i and the B axis in a regu- 
lar fashion, resulting in the aforementioned chessboard 
pattern. The V g i period (~10 mV) is set by the energy 
required for adding an extra electron to the dot (addition 
energy), whereas the B period (~0.1 T) corresponds to 
adding a flux quantum to the effective dot area. Based 
on the temperature dependence of G in the high valley 
conductance regions (not shown here), we can ascribe the 
enhancement of G to the Kondo effect. The transition 
from low to high valley conductance is associated with an 
abrupt jump of the V g \ position of the Coulomb peaks. 
In some cases (see ellipses in Fig. lb) the jump is accom- 
panied by a suppression of the peak height. Below we 
present a quantitative description for the experimental 
features discussed above, using first an intuitive, though 
quantitative, approach followed by a fully self-consistent 
simulation of our device. 

In the region where the two lowest LLs, labeled LLO 
and LL1, are occupied, the outline of the chessboard pat- 
tern can be considered to be a DD (Fig. 2a) stability di- 
agram, where now B acts as one of the gates. B couples 
to both "dots" (i.e. LLs) but with different "lever arms." 
Similarly, V g i couples with different capacitances to the 
two LLs. We write an elementary CI functional for the 



total energy, E, 

nm a 

^j[Cu(N e - C g0 V gl f + C 00 (N ie - C gl V gl ) 2 } 
2D 

- jy(N e - C g0 V gl )(N ie - C gl V gl ), 

where Afyi) is the number of electrons in LL0(1), 
e nrn = (2n + \m\ + 1)Huj + mhw c are the Fock-Darwin 
(FD) energy levels [21], and the sum is over the lowest 
two LLs. A LL consists of all n, m that satisfy LL in- 
dex A = n + | (|^| + to) = constant (i.e. or 1 for our 
specific case), and a is the spin index (since at the mag- 
netic fields considered here the Zeeman splitting is much 
smaller than the other relevant energy scales, we assume 
spin-degenerate states). Here, Co 2 = coq + uo 2 with uoq the 
bare confining frequency and co c the cyclotron frequency. 
Cij are the capacitance matrix elements (off-diagonal el- 
ements are always negative [17,22]), D = CooCn — C^. 
Here V g \ is defined relative to the gate voltage that in- 
duces ^Vq(i) e l ec t rons (i-e. ^o(i) ^ s tne num ber of electrons 
in LL0(1) at V gi =0). 
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FIG. 2. (a) Schematic of the dot in terms of two Landau 
levels (LLs). (b) Results of numerical minimization of Eq. 1. 
Color stripes are regions of constant N = No + Ni. Num- 
bers in parentheses show No , Ni . Black lines bound regions 
of constant No. Dotted lines show U N± boundaries". The 
Ni boundary lengths alternate (compare to Fig. lb). Capaci- 
tances and dot parabolicity estimated for iV~o = 25, N® — 10, 
Bo = 0.6 T: Coo = 3.3, Cu = 2.2, C Q1 = -2.1, C g0 = -0.10 
and C g i = —0.056 (all in aF), c^o = OAmeV. The dashed line 
is the approximate 3rd LL filling boundary, (c) Schematic 
honeycomb structure in the B,V g i plane resulting from the 
double dot model. Typical occupancy configurations (energy 
vs. position) of LL1 (gray arrows) and LLO (black arrows) at 
the dot boundary are shown for different honeycomb cells. 
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First, we estimate capacitances Cij and bare confining 
potential uq from the full SDF calculations for (N§, N®) 
electrons and for a "central" magnetic field Bo (see cap- 
tion, Fig. 2) [23]. Then, we numerically minimize Eq. 1 
with respect to No and TVi, keeping dj and ljq constant 
(even though in principle these change slightly over our 
(B,V g i) range). We find a "honeycomb" structure char- 
acteristic for DDs [17]. As shown in Figs. 2b and c, this 
structure clearly emerges and compares reasonably well 
with the experimental data of Fig. lb. 

Figure 2c schematically shows the population of the 
lowest two LLs in four neighboring honeycomb cells. 
B principally induces depopulation of LL1 to LLO; V g \ 
changes the total electron number N = 7V + N\ . Note 
that the dotted boundaries between N\ and N\ + 1 ( "7V~i 
borders"), are considerably shorter than the solid bound- 
aries between No and No + 1. Weak tunnel coupling of 
the inner LL (LL1) to the leads causes these Ni borders 
to appear as gaps (or small offsets) in the Coulomb os- 
cillations [24], which are clearly seen in Fig. lb. Analysis 
of Eq. 1 shows that the length of the Ni borders, in both 
B and V g i, is proportional to C\\ — |Cbi|. Therefore, the 
short N\ borders indicate that LLO strongly screens (via 
Cbi) the inner LL1. In other words, most of the self- 
capacitance of LL1 is taken up by capacitance to LLO. 
This screening of the inner LL is ultimately what pro- 
duces the chessboard pattern of the Kondo conductance 
(see below). Furthermore, the length of the N\ borders 
in Figs. 2b, c alternates, as observed in the experimen- 
tal data (Fig. lb). The shorter N\ borders appear as 
only small discontinuities in the experimental Coulomb 
oscillations. From Eq. 1 it follows that this alternat- 
ing pattern results from the spin degeneracy of the FD 
levels. Specifically, increasing B and transferring a LL1 
electron to LLO, at fixed TV, costs an additional, LLO 
level spacing when No changes from even to odd. The 
main feature of the experimental data that is not present 
in Figs. 2b, c is the increased valley conductance due to 
the Kondo effect. We show below, using our SDF analy- 
sis, that this increased valley conductance appears in the 
hexagons where No is odd. 

Our SDF calculations for realistic, 3D, lateral semi- 
conductor quantum dot structures have been thoroughly 
described before [25]. We compute self-consistent eigen- 
values £ pai eigenfunctions ^ pcr , occupancies n pa and tun- 
neling coefficients 7 pcr , with p the orbital and a the spin 
indices, as well as the total interacting energy F of the 
dot-gate- leads system, all as a function of TV, V g i and 
B. However, we cannot compute the coherent Kondo- 
assisted conductance of the dot from the ground state 
properties provided by the SDF calculation. Instead, we 
introduce an effective Kondo parameter, which enables 
us to reproduce the conductance modulation character- 
istic for the chessboard pattern. 

In Ref. [12] it was qualitatively argued that the alter- 
nating Kondo conductance with varying B at fixed N 



observed in their experiments in the strong edge state 
regime resulted from Coulomb regulated redistribution 
of electrons, one at a time, from LL1 to LLO. Since LLO 
was assumed much more strongly coupled to the leads, 
the Kondo effect was argued to occur only when No is 
odd. Here, the computed electronic structure allows us to 
quantitatively exhibit this depopulation process, to show 
how it varies in phase from one Coulomb valley to the 
next and to show how the Independent tunneling coeffi- 
cients affect the structure of individual chessboard fields. 
The model's assumption of spin-degenerate states is ac- 
counted for at the end of the paper. 

In the standard derivation of the s-d model from the 
Anderson model [26], a Schrieffer- Wolff transformation 
is used to define an effective coupling between spin- 
degenerate states of the dot (or impurity) via virtual elec- 
tron exchange with a neighboring Fermi sea. These co- 
tunneling processes then lead to coherent hybridization 
of the degenerate states, which results in a logarithmic 
increase of the conductance at low temperature, a fea- 
ture of the Kondo effect in quantum dots. In general, a 
dot can contain many singly occupied electron states, for 
which spin-flip can occur via co-tunneling contributing to 
the Kondo coupling. In order to estimate the contribu- 
tion of all singly-occupied states to the Kondo effect, we 
define a parameter K, which we call Kondo parameter, 
as the sum of all co-tunneling amplitudes that leave the 
ground state of the system unchanged except for the flip 
of a single spin 

K — ^ ^ Tip, a (1 Tip, a) \J^ip,a^ip,a 
p,a 



F(N+l,V g iB)-F(N,V g i,B) 



F(N-l,V gl ,B)-F(N,V gl ,B) s 

where a is the spin opposite to a. We conclude from the 
experimental data that the spin polarization of the leads 
is negligible, in contrast to Ref. [24]. 

A color scale plot of K in the (B,V g i) plane is shown 
in Fig. 3, for N= 32 to 39. Here, the V g i dependence 
has been approximated as follows. The two denomina- 
tors in Eq. 2 are additions energies, which are, to a good 
approximation, linear in V g \ and vanish at the Coulomb 
oscillations. Thus, F(N + l,V gU B) - F{N,V gU B) « 
(e 2 /2C) + ea(V g i - V g f min ) + e N+1 , where (e 2 /2C) = 
(d 2 F/dN 2 ) defines C, and dF(N,V g f min ,B)/dN = 
defines the valley center gate voltage, y^ m%n ? anc [ where 
& = C dot- gate /C is the lever arm (a similar analysis holds 
for the second energy denominator). By calculating the 
full electronic structure only near the valley centers, we 
can determine the CI parameters V^' mm , C and a and 
thereby show the V g \ dependence of K due to the 
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FIG. 3. Color scale plot of Kondo parameter K in the B,V g i 
plane. Red (blue) corresponds to large (small) K. Third LL 
filling begins below ~ 0.5 T. K is calculated in the Coulomb 
valley center at each B (see text). Upper panel shows Fermi 
energy vs. B for N — 37 and fixed V g i = -407 mV. Ef drops 
at each reconstruction (i.e. at the end of each Kondo zone). 
Inset, K fully calculated from Eq. 2 on fine mesh for small 
B,V g i region (i.e. no approximation for the V g i dependence is 
used here). 



addition energies, Fig. 3. The full calculation of Eq. 2 on 
a mesh of B,V g i values, which is numerically taxing, is 
shown for a small region in the inset to Fig. 3. 

The alternating pattern of Kondo "zones" in successive 
Coulomb valleys is evident. Even in this low magnetic 
field regime, the coupling of the LL0 states to the leads 
typically exceeds that of the LL1 states by two orders of 
magnitude. Therefore, although the parameter K is a 
sum of all possible co-tunneling amplitudes, the ampli- 
tudes of the LL0 states dominate. Hence, even when N 
is odd, K is negligible as long as No is even. 

Within a Kondo zone an abrupt increase of K is fol- 
lowed by a gradual decrease. This results from the con- 
traction, with £?, of the half-filled orbit at the dot edge, 
and the resulting decrease of its tunnel coefficient. When 
another electron depopulates from LL1 to that orbit, that 
spin flip process is no longer available, and K collapses. 
Depopulation of LL1 coincides with a drop of the dot's 
Fermi level [28] relative to the leads (Fig. 3, top panel). 
The denominators in Eq. 2 are responsible for the in- 
crease of K away from the valley centers; a feature which 
is clearly observed in the experiment Fig. lb. 

The approximation used in Fig. 3 includes the V g i de- 
pendence of the charging energies, but not that of the 
n Py(T . The latter dependence, specifically that of the 
transfer from LL1 to LL0, controls the slant of the chess- 
board fields. This slant only emerges in the calculation 
of K without approximation, as shown for a small region 
of V g i and B in the inset to Fig. 3. 

Finally, we comment on our assumption of spin- 
degenerate states. For N odd, this only requires that the 
Zeeman splitting (~ 10 fieV) is negligible, which is the 
case in the B range considered. However, for N even, a 



singly filled orbital in LL0 implies another one in LL1 and 
these electrons have, in general, an exchange interaction. 
In Ref. [27] it was found that this exchange interaction 
diminishes rapidly with TV. For our dot, we explicitly de- 
termine the singlet-triplet splitting for various (even) N 
and B by computing separate ground states constrained 
to S = and 5 = 1. We find that the splitting is typi- 
cally tens of fieV. The experimental signature of a split 
ground state is a split Kondo resonance (in source-drain 
voltage). Our analysis suggests that this splitting, which 
has been observed [9,15], would be characteristic of N 
even in regimes where Zeeman energy is small. 

We acknowledge financial support from the DARPA 
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